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PREFACE 
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V31S-04A,  and  V31S-A3A.  The  authors  of  this  report  were  G.  D.  Amey,  Jr.,  and  C.  W. 
Pender,  Jr.,  ARO,  Inc.  The  manuscript  (ARO  Control  No.  ARO-VKF-TR-76-1 1 1)  was 
submitted  for  publication  on  September  29,  1976. 


1 


AEDC-TR-76-148 


CONTENTS 


. 1.0  INTRODUCTION  7 

2.0  REVIEW  OF  THE  CONCEPT 7 

3.0  POSITIONAL  FORCE  EXPERIMENTS 

3.1  Apparatus  9 

3.2  Procedure  10 

3.3  Results  and  Discussion  11 

4.0  DAMPING  FORCE  EXPERIMENTS 

4.1  Apparatus  14 

4.2  Procedure  15 

4.3  Results  and  Discussion  15 

5.0  CONCLUDING  REMARKS  18 

REFERENCES  19 


ILLUSTRATIONS 


Figure 


1.  System  Concept  21 

2.  Null- Flux  Field  22 

3.  Flight  Nomenclature 23 

4.  Sketch  of  Positional  Force  Apparatus 24 

5.  Ip  Coil  Geometry  25 

6.  Photograph  of  Positional  Force  Apparatus  26 

7.  Scheme  of  Positional  Force  Apparatus  27 

8.  Photograph  of  Conductive  Shells  for  Positional 

Forces  28 

9.  Magnetic  Field  of  Line  Currents  versus 

Field  of  Current  Segments  29 

10.  Repulsive  Forces  versus  Theory 

a.  1,750  rpm  30 

b.  1,500  rpm  31 

c.  1,250  rpm  32 

d.  1,000  rpm  33 

e.  750  rpm  34 

f.  500  rpm  35 

g:  350  rpm  36 

h.  250  rpm  37 


3 


A EDC-TH  -76-1 48 


Figure  Page_ 


1 1 . Ratio  of  Repulsive  to  Drag  Forces 

a.  r/a  = 1.55  38 

b.  r/a  = 1.92  39 

c.  r/a  = 2.53  40 

12.  Ratio  of  Repulsive  to  Drag  Forces 

a.  h/a  = 6 4^ 

b.  h/a  = 4 42 

c.  h/a  = 2 43 

13.  Magnetic  Moment  Arm  versus  Con-elation 

Parameter  44 

14.  Effect  of  Misalignment  on  Positional  and 

Drag  Forces  45 

15.  Effect  of  Misalignment  on  Magnetic  Moment  Arm  46 

1 6.  Sketch  of  Damping  Force  Apparatus  42 

17.  Photograph  of  Damping  Force  Apparatus 48 

18.  Photograph  of  Conductive  Shells  for  Damping 

Forces  49 

19.  Representative  Damping  Pattern  of  Pendulum  (Strip  Chart) 50 

20-  Radial  Damping  Coefficients 

a.  Thick  Shell,  0 - ir/8  51 

b.  Thick  Shell,  0 = ir/4  51 

c.  Thin  Shell,  0 = ir/4  ...51 

21.  Angular  Damping  Coefficient 

a.  Thick  Shell,  0 = ir/8  52 

b.  Thick  Shell,  0 = tr/4  52 

c.  Thin  Shell,  0 = ir/4 52 

22.  Combination  Damping  Coefficient 

a.  Thick  Shell,  0=0  53 

b.  Thick  Shell,  0 = ir/ 16  53 

c.  Thick  Shell,  0 = tt/8  54 

d.  Thick  Shell,  0 = 3tt/16  54 

e.  Thin  Shell,  0 = 0 55 

f.  Thin  Shell,  0 = */ 16  55 

g.  Thin  Shell,  0 = ir/8  56 

h.  Thin  Shell,  0 = 3ir/ 16  56 

23.  Effect  of  Misalignment  on  Damping  Coefficient  . 57 


4 


AEDC-TR-76-1  48 


Page 

APPENDIXES 

A.  DERIVATION  AND  EVALUATION  OF  SHORT  SEGMENT 

CORRELATION  FACTOR  59 

B.  METHODOLOGY  FOR  EXTRACTION  OF  EXPERIMENTAL 

DAMPING  COEFFICIENTS 65 

NOMENCLATURE  69 


5 


AE  DC-T  R -76-1 48 


1.0  INTRODUCTION 

The  results  of  two  sets  of  experiments  conducted  in  support  of  a 
magnetic-guideway/guided-projectile  system  concept  for  aeroballistic  ranges  are  presented. 
A description  and  the  theoretical  aspects  of  the  concept  have  previously  been  presented 
in  Ref.  1.  The  intent  of  the  experimental  program  has  been  to  simulate  the  conditions 
under  which  the  projectile  is  influenced  by  the  guideway  and  to  measure  the  simulated 
"guidance”  forces  for  verification  of  theory. 

In  one  set  of  experiments,  the  objective  was  to  simulate  the  conditions  under  which 
the  projectile  experiences  a positional  force,  i.e.,  that  component  of  the  guidance  force 
that  acts  toward  the  guideway  axis  in  an  amount  proportional  to  displacement  therefrom. 
Such  positional  forces  aie  generated  electrodynamically  by  virtue  of  the  projectile's  rapid 
downrange  velocity  past  positional  elements  of  the  guideway.  Simulation  of  the  necessary 
conditions  has  been  accomplished  to  a reasonable  degree  by  placing  scaled-down  positional 
guideway  elements  near  the  periphery  of  a large  wheel  and  rotating  the  wheel  at  high 
speed.  This  effectively  creates  a traveling  guideway  with  a stationary  projectile  so  that 
simulated  positional  forces  might  be  measured  with  fixed  equipment.  Lack  of  complete 
simulation  made  it  necessary,  however,  to  develop  and  apply  a correlation  factor  to  the 
idealized  theory  of  Ref.  1 so  that  meaningful  comparisons  could  be  made  with  the 
experimental  results.  Once  this  is  done,  theory  and  experiment  agree  quite  well. 

In  the  other  set  of  experiments,  the  objective  was  to  simulate  the  conditions  under 
which  the  projectile  experiences  lateral  damping  forces,  i.e.,  that  component  of  the  guidance 
force  that  tends  to  retard  or  damp  lateral  motion  to  achieve  stability.  In  conjunction 
with  guideway  damping  elements,  these  forces  are  generated  electrodynamically  as  a result 
of  the  lateral  motion  itself.  It  has  been  possible  to  simulate  the  conditions  with  a pendulum 
and  a short  section  of  simulated  guideway.  Damping  forces  have  been  inferred  from  the 
rate  at  which  oscillations  of  the  pendulum  decay.  The  experimental  results  thus  obtained 
agree  well  with  Ref.  1 theory. 


2.0  REVIEW  OF  THE  CONCEPT 

The  general  overall  arrangement  of  the  system  concept  is  shown  in  Fig.  1.  The 
guideway  is  provided  by  four,  range-length  line  currents  of  magnitude  Id  and  four, 
range-length  sets  of  transposed  line  current  segments  of  length  (o)  and  magnitude  (Ip); 
all  of  which  are  symmetrically  orientated  about  the  range  centerline.  Also,  the  projectile 
is  a flared  cylinder  with  a conductive  shell  around  the  cylindrical  portion.  The  nose  piece 
to  be  tested  is  optional. 
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The  pattern  of  the  guideway 's  magnetic  field  produced  by  the  ID  currents  alone 
is  shown  in  Fig.  2.  A similar  so-called  null-flux  field  is  produced  by  the  Ip  currents  alone 
except  that  the  flux  is  of  opposite  direction  in  adjoining  segments.  Because  Ip  is  greater 
than  Ip,  the  resultant  field  of  all  the  currents  is  a range-length  null-flux  field  that  is 
modulated  with  respect  to  distance  downrange.  Thanks  to  the  principle  of  superposition, 
however,  that  part  contributed  by  Id  and  that  part  contributed  by  Ip  can  be  considered 
independently  insofar  as  interaction  with  the  projectile  is  concerned.  For  consideration 
of  either,  an  important  point  that  should  be  stressed  is  that  the  field  is  strong  near  the 
currents  and  decreases  to  zero  or  the  null  value  at  their  axis  of  symmetry;  i.e.,  the  guideway 
axis  which  coincides  with  the  range  centerline. 

Now  visualize  the  projectile  to  be  in  high-speed  downrange  motion  parallel  to  the 
guideway  axis.  It  will  experience  no  change  in  the  field  of  the  ID  currents  so  that  no 
electrodynamic  interaction  with  them  will  occur.  However,  as  it  passes  from  segment  to 
segment,  there  is  exposure  at  relatively  high  frequency  to  a field  of  opposite  polarity 
insofar  as  the  Ip  currents  are  concerned.  At  each  change  in  polarity,  the  intent  is  that 
the  induced  voltages  eliminate  an  existing  set  of  eddy  currents  in  the  conductive  shell 
and  create  ones  of  opposite  direction  that  will  persist  until  the  next  change.  From  Lenz' 
law,  the  field  of  the  eddy  currents  is  in  direct  opposition  to,  and  tends  to  cancel  or 
nullify,  the  change  which  caused  them.  Hence,  the  shell  tends  to  appear  as  a diamagnetic 
object  in  the  field  of  the  Ip  currents.  To  the  extent  this  is  achieved,  the  shell  is  repulsed 
by  these  currents  and  causes  it  to  seek  a position  of  equilibrium  in  the  null-flux  region 
along  the  guideway  axis.  Because  these  forces  and  their  resultant  depend  on  the  shell's 
lateral  position  with  respect  to  this  axis,  they  are  called  magnetic  positional  forces  and 
are  denoted  by  F^  as  in  Ref.  1.  Correspondingly,  the  Ip  currents  can  be  considered  as 
the  positional  force  elements  of  the  guideway. 

Also,  consider  that  there  is  a relatively  slow  lateral  motion  of  the  projectile  as  it 
progresses  downrange.  The  timewise  average  of  the  added  electrodynamic  forces  acting 
on  the  projectile  are  of  interest  to  produce  stability.  To  determine  this  average,  it  is 
necessary  to  consider  only  the  average  magnetic  field,  which  is  the  field  of  the  ID  currents. 
As  can  be  seen  from  Fig.  2,  lateral  motion  in  this  field  results  in  the  projectile's  exposure 
to  changing  flux  so  that  again  eddy  currents  are  induced  in  the  conductive  shell.  However, 
as  stated  above,  lateral  motion  is  slow  so  that  the  eddy  current  energy  is  dissipated  in 
the  shell  as  joulean  heat  at  the  same  rate  it  is  generated.  Therefore,  the  lateral  motion 
generating  the  currents  is  retarded  or  damped.  Hence,  the  ID  currents  can  be  considered 
as  guideway  damping  elements  that  apply  a damping  force  (FMD)  to  the  projectile  (Ref. 
1). 


As  indicated  in  Fig.  3,  FM  and  FM  D effectively  act  on  the  projectile  through  a 
magnetic  force  center  (cM)  that  is  a distance  (£M)  forward  of  its  center  of  gravity  (cg). 
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Also,  for  stability  purposes,  the  aerodynamic  forces  act  through  an  effective  center  of 
pressure  (cp)  at  a distance  (£)  aft  of  cg.  At  any  time  the  projectile  is  off  course,  FM 
forces  its  nose  toward  the  guideway  axis.  A change  in  angle  of  attack  is  then  to  enlist 
aerodynamic  forces  to  aid  in  correction  of  the  course.  Once  this  corrective  action  is 
underway,  Fm  d is  effective  in  reducing  overshoot  and  damping  oscillations  about  the 
guideway  axis. 


3.0  POSITIONAL  FORCE  EXPERIMENTS 


3.1  APPARATUS 

The  apparatus  designed  and  built  for  this  set  of  experiments  is  illustrated  in  Fig. 
4.  From  a mechanical  viewpoint,  it  consists  essentially  of  a large  diameter  coil  containment 
wheel  that  has  a safety  brake  and  is  rotated  by  a variable-speed  drive  unit.  From  an  electrical 
viewpoint,  it  consists  of  an  Ip  simulation  coil,  a 0-  to  1,000-amp  variable  d-c  power  supply 
and  slip  rings  for  energization  of  the  coil,  and  models  of  projectile  conductive  shells.  As 
indicated,  the  coil  is  contained  by,  and  moves  with,  the  rim  of  the  wheel  while  the  model 
is  fixed  on  a multi-component  wind  tunnel  force  balance.  The  balance  in  turn  is  mounted 
on  a stationary  cross-feed  that  has  a vernier  for  model  positioning. 

The  geometry  of  the  coil  itself  is  shown  in  Fig.  5.  It  consists  of  five  layers  or  turns 
and  provides  two  zigzag  patterns  of  current  flow  in  opposite  directions  around  the  inside 
diameter  of  the  wheel  rim.  Alternate  segments  of  the  zigzag  patterns,  shaded  in  the  sketch, 
form  a circle  around  the  rim's  axial  center,  and  it  is  these  portions  of  the  pattern  that 
are  intended  to  simulate  Ip  segments  of  the  guideway.  The  unshaded  portions  of  the 
patterns  ideally  would  be  sufficiently  remote  from  the  model  to  be  ineffective  in  generating 
forces.  The  extent  to  which  this  has  not  been  achieved,  however  (see  Fig.  6),  is  believed 
to  be  fair  compensation  for  the  fact  that  there  is  effectively  a significant  space  between 
adjoining  segments  of  interest  (in  addition  to  insulation  space  between  segments,  there  is  a 
tendency  for  currents  to  take  shortcuts  around  comers  at  their  ends). 

The  overall  idealized  scheme,  which  the  apparatus  attempts  to  emulate,  is  illustrated 
in  Fig.  7.  Ideally,  Ip  elements  would  be  made  to  travel  in  a circle,  and  each  element 
should  exist  to  the  next  element.  The  cylindrical  shell  would  then  be  positioned  a distance 
(rj ) from  the  circle  and  have  an  axial  curvature  to  match.  With  various  rotational  speeds, 
the  electrodynamic  forces  acting  on  the  cylinder  would  then  be  measured.  These  would 
include  the  force  of  repulsion  (Fm),  the  drag  force  (Fdrg)>  and  the  moment  about  the 
point  of  measurement.  With  a circle  of  sufficient  size,  the  interaction  should  approach 
that  which  would  occur  between  the  conductive  shell  on  a projectile  and  one  co-linear  set  of 
line  current  segments.  The  results  for  four  such  sets  in  a null-flux  configuration  would  be 
obtained  by  superposition. 
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For  the  present  experiments,  which  utilize  currents  that  are  distributed  over  the  finite 
cross-sectional  area  of  five  turns,  an  equivalent  filamentary  current  (Ip)  is  assumed  to 
be  at  the  geometric  center  of  this  area.  Its  effective  magnitude  is  the  total  current  in 
the  five-series  turns;  i.e.,  Ip  = 5 Is  where  Is  is  the  power  supply  current. 

The  three  cylindrical  shell  models,  which  were  utilized  in  the  experiments,  are  shown 
in  Fig.  8.  They  are  identical  except  for  having  individual  lengths  of  0.0508,  0.1016.  and 
0. 1 524  m,  respectively.  All  are  made  of  OFHC®  copper,  have  an  outside  diameter  of  0.0508 
m,  and  a wall  thickness  of  0.00254  m except  for  the  end  closures  which  are  0.00635-m 
thick.  Their  curvature  of  axis  is  intended  to  match  the  curvature  of  the  rotating  coil. 

A printer-scanner  and  an  integrating  digital  voltmeter  were  provided  for  recording 
the  experimental  data  in  digital  form.  Items  recorded  included  the  rotational  speed  of 
the  wheel  as  sensed  with  a tachometer,  the  current  Is  as  measured  using  a current  shunt, 
the  force  balance  outputs,  and  the  excitation  voltage  to  the  strain-gage  bridges  of  the 
balance. 

3.2  PROCEDURE 

Before  measuring  the  forces  on  each  of  the  three  models,  a careful  test  stand  alignment 
of  the  shell  with  the  force  balance  components  was  made,  and  they  in  turn  were  aligned 
with  the  Earth's  vertical  and  horizontal.  The  balance  was  then  calibrated  against  weights 
of  known  value. 

With  the  model  remaining  fixed  to  the  balance,  the  two  were  subsequently  moved 
to  the  rotating  apparatus  where  care  was  taken  in  attaching  the  balance  to  the  cross-feed 
to  maintain  vertical  and  horizontal  orientation.  The  cross-feed  was  then  raised  or  lowered 
to  provide  best  alignment  between  the  model  contour  and  the  coil.  A chosen  reference 
distance  (r! ) between  the  shell  axis  and  the  coil  center  was  then  carefully  established, 
and  a corresponding  reference  reading  on  the  cross-feed  vernier  was  recorded. 

To  begin  the  simulation  and  measurement  of  forces,  a desired  rotational  speed  was 
first  established,  and  the  model  was  positioned  at  a desired  r!  through  appropriate  setting 
of  the  cross-feed  vernier.  Finally,  with  zero  conditions  being  recorded,  both  before  and 
after,  the  Ip  coil  was  energized  for  a period  of  30  sec,  during  which  force  measurements 
were  made.  The  coil  was  then  allowed  to  continue  rotating  while  cooling  for  a period 
of  five  minutes  prior  to  the  next  data  point.  The  process  was  continued  until  a grid  of 
data  was  obtained  for  speeds  from  250  to  1,750  rpm  and  for  values  of  from  0.037 
to  0.072  m.  During  the  process,  some  additional  data  were  taken  with  the  model  of  medium 
length  to  determine  the  effect  of  misalignment  between  model  and  coil. 
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Finally,  to  enable  a more  accurate  determination  of  ri  under  actual  test  conditions, 
it  was  necessary  to  measure  the  containment  wheel  growth  versus  speed  and  the  force 
balance  deflection  versus  force.  As  might  be  expected,  the  wheel  growth  was  found  to 
be  proportional  to  the  square  of  its  rotational  speed,  and  the  balance  deflection  proved 
to  be  a linear  function  of  the  applied  forces. 

3.3  RESULTS  AND  DISCUSSION 

The  Ref.  1 theory,  which  the  experimental  results  are  intended  to  validate,  predicts 
that  in  an  idealized  experiment  the  force  of  repulsion  per  unit  length  on  the  conductive 
shell  should  be 


Fr  = 


/i0Ip<F  ( S - 1 ) 

2irri  (r|  2 - a) 


0) 


where 


is  the  permeability  of  free  space 

Ip  is  the  simulated  guideway  current 

a is  the  effective  shell  radius 

rj  is  the  distance  from  lp  to  the  shell  axis 

and  where  £ is  the  ratio  of  the  magnetic  field  within  the  cylinder  to  that  which  would 
exist  at  the  same  location  if  the  cylinder  were  absent.  Reference  1 gives  its  value  as 


$ 


(2) 


where  <;  is  the  area  resistivity  of  the  shell  and  co  is  the  angular  velocity  at  which  it 
is  exposed  to  a field  of  alternate  polarity  by  the  passing  lp  elements. 

Unfortunately,  certain  compromises  had  to  be  made  in  design  of  the  apparatus  so 
that  the  actual  experiments  have  been  somewhat  less  than  ideal.  As  is  seen  from  Eq.  (1), 
it  is  desired  that  > be  small  so  that  Fr  will  be  large.  From  Eq.  (2)  then,  c must  be 
small,  or  the  product  of  to  and  a must  be  large.  In  fact,  however,  the  radius  (a)  had 
to  be  limited  to  a relatively  small  value  to  match  the  size  of  the  rotating  apparatus.  It, 

in  turn,  was  limited  in  size  by  availability  of  resources.  Because  the  choice  was  made 
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to  avoid  cryogenics,  the  effective  <:  was  limited  to  the  relatively  large  values  which  could 
be  obtained  at  room  temperature  with  OFHC  copper  and  a shell  thickness  to  match  a. 
Thus  forced  toward  large  values  for  c o and  limited  to  relatively  low  peripheral  velocities 
achievable  with  rotating  machinery,  a small  segment  length  (Q)  was  necessary  for  the 
experimental  Ip.  This  compromise  is  considerably  less  severe  in  an  actual  range  where 
projectile  velocities  are  4,000  to  6,000  m/sec. 

When  compared  with  all  radial  guideway  and  projectile  dimensions,  Ref.  1 theory 
presumes  that  a can  be  at  least  20  times  greater  and  thus  treats  segments  as  being  infinite 
in  deriving  theory.  The  implications  of  a not  being  so  great  was  made  clear  by  computing 
the  ratio  (Btrp/Bjnf)  as  a function  ofz'  for  various  values  of  a/ r’  where  Bjrp  is  the 
field  that  would  exist  a distance  (r')  from  a co-linear  set  of  transposed  segments  of  length 
(«)  and  where  Bjnf  represents  that  field  which  would  exist  if  the  set  were  replaced  by 
an  infinite  Line  current.  A plot  of  the  computed  results  (Fig.  9)  shows  that  the  average 
field  intensity,  to  which  a conductive  shell  is  exposed,  is  weak  for  o/r'  ^ 5 as  was  the 
case  in  the  present  experimental  situation.  It  is,  therefore,  necessary  to  develop  and  apply 
a correlation  factor  (fCon(a,  rj ))  to  the  Ret.  1 theory  for  meaningful  comparisons  to 
be  made  with  the  experimental  results.  The  derivation  of  expressions  for  this  factor  is 
given  in  Appendix  A,  which  includes  a computer  program  for  evaluation  under  specific 
conditions. 

The  actual  experimental  results  for  Fr  together  with  Ref.  1 theory  and  Ref.  1 theory 
times  fc o r r , a,  rx ),  as  the  latter  applies  to  the  experimental  situation,  is  shown  in  Figs. 
10a  through  h.  The  data  show  that  there  is  no  consistently  discernible  difference  in  the 
results  for  the  three  different  h/a  ratios  tested.  This  is  as  it  should  be  since  Fr  is  force 
per  unit  length.  However,  with  thick  end  closures,  the  h/a  ratio  is  not  required  to  be 
excessively  great  for  two-dimensional  mathematical  modeling  to  be  satisfactory. 

The  data  are  considerably  below  predictions  of  the  idealized  Ref.  1 theory,  especially 
at  increasing  values  of  ri/a.  The  percentage  of  disagreement  is  considerably  less  at  small 
fi/a,  as  was  expected.  The  application  of  fcorr  (a,  a,  rx)  to  the  theory  brings  it  into 
harmony  with  the  data  at  the  larger  values  of  rx  / a,  but  overcompensates  at  the  smaller 
values  of  this  ratio.  One  explanation  of  this  is  that  the  effective  location  of  1P  is  closer 
to  the  shell  than  was  assumed;  i-e.,  the  effective  rt  may  be  less.  If  it  is  considered  that 
geometric  mean  distances  are  normally  used  to  quantitate  the  mutual  inductance  between 
conductors,  this  explanation  is  quite  plausible  (Ref.  2). 

Another  possible  reason  for  the  discrepancy  is  that  the  effective  value  of  a = 0.02413 
m which  was  used  to  correlate  the  data  was  not  the  proper  choice.  This  was  based  on 
the  distance  from  the  shell  axis  to  the  midpoint  of  its  wall.  Use  of  a = 0.0254  m (the 
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radius  of  the  outside  surface)  would  have  resulted  in  better  agreement.  Neither  explanation 
will  be  pursued,  however,  as  the  agreement  shown  is  good  considering  the  compromises 
in  both  theoretical  derivations  and  experimental  design.  Agreement  is  adequate  for 
rotational  speeds  of  the  apparatus  from  250  to  1,750  rpm  where  theoretical  values  of 
6 from  0.67  to  0.13  were  attained. 

In  addition  to  confidence  in  predicting  values  for  Fr,  it  is  important  that  its  value 
be  large  with  respect  to  the  corresponding  drag  force  (Fdrg).  Therefore,  experimentally 
measured  values  of  F^rg  have  been  utilized  in  conjunction  with  the  experimental  values 
of  Fr  to  formulate  the  ratio  (Ft/FDRG),  The  results  for  the  three  shell  models  at  four 
different  n/a  positions  are  plotted  in  Figs.  I la  through  c and  12a  through  c,  inclusively. 
From  these,  it  is  clear  that,  as  the  correlation  parameter  (^0  a co/ 2c)  (see  Eq.  (2))  is 
made  small.  Fdrg  becomes  a significant  factor,  whereas  it  becomes  insignificant  as 
compared  with  Fr  when  this  parameter  is  made  large.  This  trend  was  predicted  in  Ref. 

1 based  on  the  fact  that,  for  low  values  of  the  parameter,  the  eddy  currents  in  the  shell 
are  dissipated  at  the  rate  they  are  induced,  whereas,  for  high  values  they  tend  to  be 
persistent  and  result  in  shielding.  From  Figs.  11a  through  c alone,  it  is  also  clear  that 
the  magnitude  of  Fdrg  slightly  increases  with  decreasing  values  of  the  shell  s h/a  ratio. 
While  not  serious  for  anticipated  values  of  h/a,  this  trend  is  to  be  expected  from  the 
increasing  importance  of  end  effects.  Similarly  from  Figs.  12a  through  c alone,  it  is  clear 
that  the  magnitude  of  FDrc  also  increases  slightly  with  increasing  values  of  the  rj/a 
ratio.  In  this  case,  the  trend  can  be  attributed  to  the  fact  that,  at  larger  ri  /a,  the  inducing 
field  and  the  eddy  current  field  become  less  two-dimensional  so  that  the  eddy  current 
pattern  is  more  circular.  If  * were  larger  as  would  be  the  case  in  an  actual  guideway, 
tills  trend  would  not  be  expected. 

Equally  important  to  the  magnitudes  of  the  forces  is  the  establishment  of  a point 
on  the  shell  through  which  they  can  be  assumed  to  act.  To  resolve  this,  the  moment 
was  measured  each  time,  about  the  same  point  at  which  Fr  and  Fdrg  were  being  sensed. 
For  convenience,  the  authors  have  chosen  to  use  this  measured  moment  to  transfer  Fr  and 
Fdrg  to  effective  point  of  action  that  lies  on  the  shell  axis  at  a distance  (AUm  ) aft  of 
the  shell's  geometric  center.  In  this  way,  although  the  moment  is  essentially  a result  of 
Edrc  being  off  axis,  its  effect  can  be  accounted  for  by  a change  (A£M)  in  a guided 

projectile's  magnetic  moment  ami  (JEm  ).  It  so  happens  that  the  experimental  results  collapse 

« 

to  a simple  function  of  the  correlation  parameter  (ja0  a co/2  O as  shown  in  Fig.  13.  At  larger 
values  of  this  parameter  where  a guideway  system  would  operate,  the  magnitude  of  (Afi^ ) 
decreases. 
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A final  set  of  positional  force  data  was  taken  to  determine  the  effects  of  misalignment 
between  the  shell  axis  and  Ip,  as  might  be  expected  in  the  case  of  a projectile  in  a guide  way. 
This  was  done  at  a relatively  large  fixed  value  of  the  correlation  parameter  and  at  a fixed 
rt/a  ratio  with  the  angle  included  between  Ip  and  the  shell  axis  being  varied  up  to  10 
deg.  The  results  for  Fr  and  FDRC  versus  the  angle  are  shown  in  Fig.  14,  where  it  can 
be  seen  that  they  are  little  perturbed.  However,  there  is  a significant  effect  insofar  as 
their  effective  location  is  concerned  as  illustrated  in  Fig.  15  where  A6M/a  is  plotted  versus 
the  angle.  Care  must,  therefore,  be  taken  when  designing  actual  projectiles  for  guidance 
to  ensure  that  £M  /a  is  sufficiently  large  and  that  the  included  angle  between  a projectile 
axis  and  the  guideway  axis  does  not  assume  excessive  values  during  flight. 

4.0  DAMPING  EXPERIMENTS 


4.1  APPARATUS 

Since  electromagnetic  damping  of  a projectile's  lateral  motion  is  a consequence  of 
only  that  component  of  motion  which  is  normal  to  the  guideway  axis,  it  was  possible 
to  devise  an  experimental  apparatus  closely  simulating  the  necessary  conditions  of  an  actual 
system.  The  apparatus  consisted  essentially  of  a pendulum  whereby  a cylindrical  copper 
shell  or  model  with  closed  ends  resembling  the  conductive  portion  of  the  guided  projectile 
in  Ref.  1 was  suspended  in  a null-flux  field  of  four  simulated  Id  currents  (see  Figs.  1 
and  2).  Figure  16  shows  a cross  section  of  the  apparatus,  and  Fig.  17  is  a photograph 
of  the  actual  hardware. 

The  lD  currents  were  simulated  by  a coil  of  nine  turns  of  copper  tubing  wound 
on  a wooden  form.  Since  the  power  supply  provided  a maximum  of  only  1,000  amp, 
a nine-turn  arrangement  was  used  to  increase  the  effective  ID  and  the  resulting  null-flux 
field  for  greater  damping,  i.e.,  ID  = 9 ls.  To  reduce  possible  end  effects  of  the  coil,  its 
length  was  made  significantly  longer  than  its  radial  dimension  (b);  a 1.2-m  length  was 
used  with  the  effective  b dimension  of  0.0898  m.  Since  a maximum  energy  dissipation 
rate  of  approximately  65,000  w was  required,  water  cooling  was  employed. 

Four  different  shell  models  (Fig.  18)  were  used  to  investigate  the  effects  of  different 
parameters.  All  had  an  outside  diameter  of  0.0254  m;  two  were  0.005  m and  two  were 
0.01  m long.  One  of  each  length  had  a wall  thickness  of  6.35  x 10’4  m,  and  one  of 
each  length  had  a wall  thickness  of  1.245  x 10’3  m.  All  had  4.775  x 10  3-m-thick  copper 
end  closures  attached  with  silver  solder. 

The  cylindrical  shells  were  drilled  through  at  their  midpoints  and  were  attached  to 
threaded  stainless  steel  rods  with  brass  nuts  (in  an  attempt  to  have  no  ferromagnetic 
materials  involved).  These  assemblies  then  had  razor  blades  secured  to  their  other  ends 
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(see  Fig,  16)  to  provide  low  friction  pivots  when  placed  in  the  grooved  fixture  on  the 
stand  holding  the  Id  coil. 

The  apparatus  was  designed  to  allow  adjustment  of  the  model's  position  in  the 
null-flux  Field  and,  thereby,  acquire  damping* data  at  various  locations. 

Data  were  taken  using  a prototype  laser  vibrometer  (Ref.  3).  Use  of  this  instrument 
allowed  continual  non-interference  monitoring  of  the  model's  position  and  velocity.  The 
models  were  equipped  with  retroreflectors  (comer  reflectors),  and  the  instrument  was 
located  approximately  one  meter  from  the  pendulum.  Data  output  was  in  the  form  of 
strip  charts  (see  Fig.  19  for  a typical  example). 

4.2  PROCEDURE 

The  procedure  was  to  place  each  model  at  some  selected  (and  known)  location  in 
the  field,  then  to  displace  it  slightly,  release  it,  and  observe  its  decaying  oscillatory  motion. 
Adjustments  on  the  stand  were  used  to  shift  the  pendulum  either  horizontally  or  vertically. 
A piece  of  clear  cellophane  with  grid  markings  was  mounted  near  one  end  of  the  pendulum 
to  facilitate  locating  the  models  at  the  desired  positions.  Data  were  then  taken  for  the 
various  models  at  numerous  positions  using  different  current  levels.  Data  were  also  taken 
with  no  current  to  determine  damping  due  to  factors  other  than  the  magnetic  field. 

Current  level  in  the  coil  was  determined  from  its  terminal  voltage  after  having  used 
a Kelvin  bridge  to  determine  its  resistance.  Coil  temperature  was  monitored  with 
thermocouples  and  a digital  meter. 

After  experiments  had  been  completed  with  the  coil  situated  as  shown  in  Fig.  16, 
the  coil  was  rotated  approximately  22.5  deg  about  its  indicated  pivot,  and  data  were 
taken  at  numerous  locations  in  the  rotated  field.  This  was  done  to  enable  measurement 
of  purely  radial  and  purely  angular  damping  coefficients  for  an  angular  position,  or  0, 
other  than  3ir/4  and  tt/4,  respectively. 

4.3  RESULTS  AND  DISCUSSION 

In  the  experimental  setup,  the  various  dimensions  were  of  the  same  order  of  magnitude 
as  in  the  proposed  system  (Ref.  1).  However,  there  were  some  quantities  which  differed 
appreciably  from  what  would  be  expected  in  an  actual  guidance  system.  The  Id  current 
levels,  despite  the  effect  of  having  a nine-tum  arrangement,  were  significantly  lower  than 
the  100,000  amp  used  in  Ref.  1.  The  lateral  velocity  of  the  models  in  the  experiments 
was  also  significantly  different  from  that  anticipated  during  guided  flight.  These  differences 
amounted  to  approximately  two  orders  of  magnitude  in  the  case  of  the  current  level  and 
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several  orders  in  the  case  of  lateral  velocity.  However,  damping  forces  were  generated 
which  were  adequate  for  measurement  and  correlation. 

The  intent  of  the  experiments  was  primarly  to  verify  magnitudes  of  the  damping 
coefficient  (K^  d)  in  the  central  region  of  the  field  where  K/b  < 0.3,  and  to  verify  trends 
in  that  coefficient  in  the  region  where  IR/b  > 0.3.  Because  theory  predicted  that  the 
damping  was  nearly  uniform  in  the  region  where  !R/b  ^ 0.3,  the  guidance  system  relied 
on  restricting  flight  to  that  region  and  using  a constant  value  for  KMD.  It  was,  therefore, 
important  that  this  constant  be  confirmed.  It  is  also  interesting  to  observe  how  well  theory 
and  experiment  agreed  in  the  outer  region  where  £R/b  > 0.3.  In  addition  to  the  inaccuracies 
of  the  simulation  already  emphasized,  filamentary  currents  were  assumed  in  the  theoretical 
treatment,  whereas  nine  turns  of  0.008-m-diam  copper  tubing  were  used  in  the  experiment. 
There  was  the  possibility  that  the  assumed  b was  not  accurate  since  the  copper  tubing 
was  somewhat  difficult  to  position  on  the  form  and  rather  bulky  insulating  material  was 
used  between  the  conductors.  These  experimental  inaccuracies  were  presumed  to  have  a 
more  dramatic  effect  as  S/b  increased,  and  therefore,  some  error  was  to  be  expected 
in  the  outer  region. 

Experimental  results  took  the  form  of  strip  charts,  (see  Fig.  19).  Data  retrieved  from 
these  charts  consisted  of  the  magnitudes  of  the  quantities  identified  as  2Xj , 2X2 , and 
At.  As  explained  in  Appendix  B (Eq.  (B-5)),  these  are  adequate  for  extraction  of  K^d 
from  the  relationship 


Kmd  = 

*M 


^ *n  >11  - CD 
At  X2 


(3) 


where 


and 


CD  = Damping  due  to  factors  other  than  magnetic  field  (Eq.(B-6)) 

M = Mass  of  pendulum 

£ = Separation  between  the  pendulum's  pivot  point  and  the  center  of  gravity. 

= Separation  between  pivot  point  and  center  of  magnetic  forces 


According  to  Ref.  1 theory,  it  should  then  be  possible  to  present  the  results  in 
normalized  form  by  plotting  K^ob  <:/(hI2)  versus  S/b  for  each  0.  For  this  reason, 
numerous  measurements  were  performed  prior  to  the  actual  tests.  A summary  of  the 
measurements  is  tabulated  below  for  use  as  needed.  Other  measured  quantities  not  shown 
include  a chart  speed  of  2.52  x 10  3 m/sec  and  a coil  resistance  of  6.58  x 10'2  ohm. 
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Is  currents  from  400  to  700  amp  were  used.  This  provided,  since  nine  turns  were  used, 
a magnetic  field  equivalent  to  one  produced  by  a filament  carrying  from  3,600  to  6,300 
amp. 


No. 

Shell 

Length, 

h,  j 

m x 1(T 

Shell 

Wall, 

w», 

m x 10^ 

Shell 

Radius, 

a”.  2 

an  ID4 

Area 

Resistivity , 

<***•  s 
ohn  x 10 

Coll 
Bore, 
b,  2 
■ x 10 

Pendulum 

length, 

l. 

m x 10 

Pendulum 

length, 

» Xo 

Pendulum 
Mass, 
M,  2 
K*  x 10* 

1 

5.0 

6.35 

1.24 

6.563 

0.98 

2.643 

3.088 

8.84 

2 

5.0 

12.446 

1.21 

2.427 

2.719 

3.119 

10.49 

3 

10.0 

6.35 

1.24 

6.563 

2.635 

3.033 

8.69 

4 

10.0 

12.446 

1.21 

2.427 

2.762 

3.119 

11.96 

♦ Side  wall  only;  end  walls  were  4.775  x 1(T3  m for  all  models 
**  Radius  based  on  wall  midpoint 

***  Based  on  measured  end-to-end  resistance  of  long  models  which  was  3.25  x 10'-'  and  8.57  x 10"3 
ohm  for  thick  and  thin  shell,  respectively. 


During  the  experiments,  it  was  found  that  the  stainless  steel  pendulum  rod  appeared 
to  be  slightly  ferromagnetic.  This  caused  some  displacement  of  the  pendulum's  equilibrium 
position  in  part  of  the  experiments.  The  displacement  was  determined  from  the  strip  charts, 
and  corrections  for  equilibrium  position  were  included  in  the  results  to  be  presented. 

Theory  developed  in  Ref.  1 indicated  that  the  damping  coefficients  arising  from  radial 
motion  or  from  angular  motion  should  be  identical.  It  was  convenient  to  simulate  purely 
radial  and  purely  angular  motion  relative  to  the  axis  of  the  null-flux  field  at  only  certain 
angular  positions  because  of  interference  between  the  coil  and  pendulum.  Note  in  Fig. 
16  that,  at  points  where  the  pendulum's  motion  is  at  right  angles  to  a radius  from  the 
coil  axis  through  the  equilibrium  position  of  the  pendulum,  the  motion  is  angular,  while 
that  motion  on  a path  coinciding  with  a radius  is  radial.  Figures  20a  through  c present 
experimental  results  and  theoretical  predictions  for  damping  coefficients  resulting  from 
the  damping  of  radial  motion,  whereas  Figs.  21a  through  c depict  results  arising  from 
the  damping  of  angular  motion.  The  results  shown  in  Fig.  20a  were  obtained  with  the 
thick  shell  model  with  |3  (see  Fig.  16)  equal  to  jt/8  and  with  $ varying  from  0 to  0.5 
b.  As  the  captions  indicate,  Figs.  20b  and  c and  21a  through  c are  for  a variety  of  situations 
with  different  shells,  locations,  and  directions  of  motion.  Examination  of  these  figures 
reveals  satisfactory  agreement  between  experimental  and  theoretical  results  and  confirms 
that  the  damping  coefficient  is  insensitive  to  whether  the  direction  of  motion  is  radial 
or  angular. 
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Since  theory  and  experiment  suggested  that  the  damping  coefficient  is  independent 
of  direction  of  motion  and  since  the  experimental  apparatus  was  more  conducive  to 
obtaining  results  at  paints  where  motion  of  the  pendulum  was  not  purely  radial  or  angular, 
numerous  results  were  obtained  for  situations  where  motion  was  a combination  of  radial 
and  angular.  Figures  22a  through  h show  the  resulting  combination  damping  coefficients 
for  the  different  shell  models  and  positions.  These  figures  include  the  theoretical  predictions 
for  damping  coefficients  obtained,  assuming  either  radial  or  angular  motion,  and  provide 
confirmation  that  direction  of  motion  is  not  an  influencing  factor  on  damping.  Satisfactory 
agreement  is  also  revealed  between  theory  and  experiment  where  K/b  is  small  and  trends 
are  shown  to  be  as  expected  for  the  larger  values. 

It  is  noticeable  that  experimental  damping  coefficients  for  the  thin  wall  models  are 
consistently  larger  than  those  for  the  thicker  wall  models.  This  is  reasonable  since  the 
theory  does  not  take  into  account  eddy  current  losses  in  the  rather  thick  end  plates. 
It  is  to  be  expected  that  these  losses  would  be  a significantly  greater  percentage  of  the 
total  in  the  thin  models  as  compared  to  the  thick  ones.  The  influence  of  the  end  plates 
is  more  evident  in  the  short  model  data  which  are  included  in  Figs.  20b,  20c,  21b,  and 
21c.  It  is  in  the  case  of  the  thin  short  models  that  the  effect  is  most  pronounced. 

Reference  1 assumed  that,  since  any  deviation  of  the  model's  axis  from  parallel  to 
the  guideway  axis  would  be  small,  the  effect  on  damping  coefficients  would  be  minor 
and  could  be  neglected.  This  assumption  was  experimentally  verified  as  shown  in  Fig. 
23.  The  results  are  presented  as  ratios  of  misaligned  versus  aligned  damping  coefficients. 
Shell  models  were  misaligned  as  much  as  14  deg  with  no  noticeable  effect  on  damping 
coefficients. 


5,0  CONCLUDING  REMARKS 

As  shown  in  Ref.  1,  the  determination  of  the  feasibility  of  an 
electromagnetic-guideway/guided-projectile  system  was  dependent  on  the  use  of 
theoretically  derived  magnetic  force  coefficients.  The  purpose  of  the  experimental  program 
described  herein  was  to  validate  those  coefficients.  Two  sets  of  experiments  were  used 
to  simulate  positional  and  damping  forces  as  would  be  generated  in  the  actual  guidance 
system.  These  were,  in  turn,  used  to  find  force  coefficients  for  comparison  with  those 
predicted  using  techniques  developed  in  Ref.  1. 

The  experimental  apparatus  used  for  generating  positional  forces  had  much  shorter 
"positional"  current  segments  than  proposed  for  the  actual  system.  This  necessitated 
developing  and  utilizing  a correlative  factor  when  comparing  the  theoretical  predictions 
from  Ref.  1 with  experimental  results.  The  comparison  clearly  indicated  the  validity  of 
the  analytical  techniques,  substantiating  both  trends  and  magnitudes. 
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In  the  case  of  the  damping  experiments,  the  apparatus  closely  resembled  the  proposed 
system  and  no  corrections  were  indicated.  Results  were  in  close  agreement  with  theory 
and  very  satisfactorily  corroborated  the  accuracy  of  the  previously  developed  theory. 
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Figure  4.  Sketch  of  positional  force  apparatus. 
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Figure  6.  Photograph  of  positional  fores  apparatus 


Figure  7.  Scheme  of  positional  force  apparatus. 
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Figure  8.  Photograph  of  conductive  shells  for  positional  forces. 
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a 1,750  rpm 

Figure  10.  Repulsive  forces  versus  theory. 
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b.  1,500  rpm 
Figure  10.  Continued. 


31 


N/funp 


AEDC-TR-76-148 


1.5  2.0  2.5  3.0 

r]/  a 

d.  1,000  rpm 
Figure  10.  Continued. 
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e.  750  rpm 
Figure  10.  Continued. 
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Figure  10.  Continued. 
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Figure  10.  Continued. 
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Figure  10.  Concluded. 
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a.  ri/a  — 1.55 

Figure  11.  Ratio  of  repulsive  to  drag  forces. 
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b.  ^ /a  = 1.92 
Figure  11.  Continued. 
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c.  r,  /a  = 2.53 
Figure  11.  Concluded. 
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Figure  12.  Ratio  of  repulsive  to  drag  forces. 
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b.  h/a  — 4 
Figure  12.  Continued. 
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c.  h/a  ^ 2 

Figure  12.  Concluded. 
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Figure  14.  Effect  of  misalignment  on  positional  and  drag  forces. 
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Figure  15.  Effect  of  misalignment  on  magnetic  moment  arm. 
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Figure  16.  Sketch  of  damping  force  apparatus. 


AEDC-TR-76-148 


Slotted  Holes  lorB 
Model  Posit  ion  irurgH 


Knife  Edge 


Pendulu 

Support 


Model 


ID  Coilf* 


Wooden  Form 


Pivot  Support 


A E D C 
5447-76 


Figure  17.  Photograph  of  damping  force  apparatus. 
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Figure  18.  Photograph  of  conductive  shells  for  damping  forces. 
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(■« — Pendulum  Displaced  and  Released 
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Figure  19.  Representative  damping  pattern  of  pendulum  (strip  chart). 
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a.  Thick  shell,  0 = jt/8 


b.  Thick  shell,  0 = tt/4 


c.  Thin  shell,  0 = n/4 
Figure  20.  Radial  damping  coefficients. 
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d.  Thick  shell,  p = 3tt/16 
Figure  22.  Continued. 
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a/b 


g.  Thin  shell.  0 = rr/8 


h.  Thin  shell,  0 = 3jt/1  6 
Figure  22.  Concluded. 
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Figure  23.  Effect  of  misalignment  on  damping  coefficient 
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APPENDIX  A 

DERIVATION  AND  EVALUATION  OF  SHORT  SEGMENT  CORRELATION  FACTOR 

Consider  a cylindrical  shell  of  radius  (a)  and  unit  length  that  has  been  moved  from 
infinity  to  a distance  (ri)  from  an  infinite  line  current  or  an  infinite  set  of  transposed 
line  current  segments  that  lie  on  a z'  axis  as  shown  in  Figs.  A-la  and  b.  Visualize  further 
that  persistent  eddy  currents  were  induced  in  the  shell  as  a result  of  this  action.  From 
Ref.  4,  the  energy  (W')  of  these  eddy  currents  in  terms  of  their  magnetic  induction  (B') 
is  given  by 

W'  - — f (B')2  dv  (A-l) 

2fl0  J » 

where  ttQ  is  the  permeability  of  free  space  and  v is  volume.  Next  let  the  cylinder  be 
moved  a distance  (drj ) with  a corresponding  change  (dW')  in  eddy  current  energy.  This 
change  in  energy  can  be  equated  to  mechanical  work  that  was  done  in  the  act,  i.e., 

dW'  = Fdrx  (A-2) 

where  F is  force  per  unit  length.  From  this,  it  is  seen  that 

F - dW'/d rt  (A-3) 


or  combining  with  Eq.  (A-l)  gives 


F = — — V (B')2  dv  (A-4) 

2Po  drj 

For  near  two-dimensional  geometry,  one  can  assume  that  the  total  field  energy  is 
proportional  to  that  within  the  shell.  Therefore  by  multiplying  by  an  assumed  constant 
of  proportionality  (K),  the  integral  of  Eq.  (A-4)  need  only  be  over  the  volume  of  the 
cylinder;  hence,  Eq.  (A-4)  becomes 

F = — -d-  f (B')2  dv  (A-5) 

2p0  drj  Jcyl 

If  there  is  to  be  a perfectly  diamagnetic  effect  (complete  shielding  inside  the  shell),  then 
for  this  space, 


or 


B' Binf 

B'  = - Bj  r p (A-6) 
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where  Bin  f and  B-jrp  are  the  fields  of  the  parallel  line  currents  or  transposed  segments 
as  the  case  may  be.  Now  let  a correlation  factor  fcorr(s  a,  rt)  be  defined  such  that 


fcort  - Ftrp/Finf 

Combination  of  Eqs.  (A- 5),  (A-6),  and  (A-7)  then  yields 


(A-7) 


icon 


(^,a  ,ri ) = < — 


l 


dri  «c  Y L 


B 


TRP 


dv 


J 


dri  J C Y L 


Bin f dv 


(A-8) 


From  Eq.  (21)  of  Ref.  1 and  with  reference  to  Fig.  (A-l), 


BTrp 


Mol 

4irr' 


z\  - n a- 


(z'j  - n a)2  + (r')2 


z\  -(n  - 1)^  ) 

4n^wt(A'9) 


where  I is  the  magnitude  of  current  in  the  line  segments,  *.  is  the  segment  length,  and 
n represents  the  nth  segment  from  the  point  where  z'  = 0.  Also  from  Eq.  (22)  of  Ref. 

1, 


Binf  “ ~ (A- 10) 

2irr 

For  any  given  z , it  is  seen  that  Binf  and  BtrP  vary  only  as  the  magnitude  of  rr  varies 
so  that  a geometry  suited  to  the  needed  integrations  is  shown  in  Fig.  A-2.  To  get  an 
average  force  per  unit  length  in  the  transposed  case,  it  is  necessary  to  average  over  one 
segment  length  (*)  at  some  relatively  large  z'  (see  also  Fig.  9).  Then  Eq.  (8)  for  the  proposed 
method  of  integration  becomes 


The  tabulation  of  a computer  program  that  has  been  used  to  solve  for  feorr(i,  a,  rt) 
is  included  in  Table  A-l.  The  program  is  truncated  at  n = 5.  A plot  of  the  computed 
results  for  a variety  of  conditions  is  shown  in  Fig.  A-3  with  those  for  the  present 
experimental  situation  where  ^/a  = 4.14  is  shown  as  the  dashed  curve.  As  shown,  fcotI(s 
a,  rx ) is  quite  significant  for  this  value  of  i/a,  especially  as  rx  /a  increases.  However,  for 
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*/a  values  on  the  order  of  sixty,  as  would  be  the  case  in  an  actual  range  guidance  system, 
fconK  a,  rj)  approaches  unity  for  expected  values  of  r^/a.  Hence,  in  most  actual  range 
applications,  the  idealized  theory  of  Ref.  1 should  be  adequate  without  fC0„U,  a,  ri)  being 
needed. 

z 


Figure  A-1.  Nomenclature  for  cylindrical  shell  and  currents. 
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Table  A-1. 


• » 
»» 

FORTRAN  PROGRAM 

«« 

«« 

•• 

FOR 

«• 

• • 
• • 

SHORT  SEGMENT  CORRELATION  FACTOR 

• • 
«« 

THIS  PROGRAM  FINDS  THE  RATIO  OF  THE  DERIVATIVES  OF  TWO  INTEGRALS 
DO  LOOP  1 EVALUATES  A SUMMATION  FOUNU  IN  THE  1NTERGRAL  IN  THE  NUMERATOR 
DO  LOOP  2 STEPS  TN  THE  l DIRECTION 
00  LOOP  3 STEPS  IN  THE  R>  OIRECTION 

DO  LOOP  4 STEPS  R 1 A S* ALL  AMOUNT  FOR  CALCULATION  OF  DERIVATIVE 

DO  LOOP  5 STFPS  TN  THE  RL  OIRECTION 

00  LOOP  100  STEPS  THF  PROGRAM  THROUGH  DIFFERENT  S TO  A RATIOS 


DIMENSION  ANUMI2) ,ADEN(J) 

C A IS  THF.  EFFECTIVE  RADIUS  OF  MODEL 
A a *024 1 3 
00  100  K-20.60.10 
AAAsK 

C S TS  THE  LENGTH  OF  A TRANSPOSED  SEGMENT 
S3AAA*A 
U«l=. 1*A 
00  5 la  U21 
AS»r-i 

Rial ,2*A»AS*DR1 
ROAaR 1 /A 
OR  1 l3.01#A 
00  4 J=1.2 
A4a  J-l 
RlaRl.A4«nRIl 
VIs.O 
V2=,  0 

0RPR=.02S«A 
RPR  = R 1 - A-DRPR 
00  3 Ma 1 *Q  1 
RPRaRPR*  DRPR 

DZ=.02S*S 

ZPRC4.0«S_DZ 
00  2 L= 1 *2 1 
ZPRaZPR.07 
FN2  = .0‘ 

00  1 N3 1*9 
A1*=N 

FN1  = ( |-l,)«»(N»ll  I *(  ( ZPR*A1»S)  /SORT  H ZPR-A1  *S  I "^RPR 
l**2> -I2PR-I A 1 — 1 . ) »S) /SORT ( (ZPR-l A 1-1. I »S1 *»2*RPR*»2) ) 
1 FN2=  FN2+FN1 


BNTRP=FN2/ (2,*RPRI 
BNINF  al ,/RPR 

C0M0N°RPR*ARC05  [ A (RPH*»2I » <fl 1 »*2 1 - | a»*2>  > / ( 2. »RPR»R1 ) ) 

Va(BNTRP*«2) »COHON«DZ*ORPR 

vl-v-vi 

Ra(BNINF**2)«C0M0N»0Z»0RPR 
2 V2=B*V2 

3 CONTINUE 
ANUMIJt-Vl 
4 ADEN  I J) aV2 

CF=  IANUM(2)-ANUM(1)J  / ( ADEN (2 ) -ADEN ( 1 ) I 

RES=S/A 

5 PRINT  20.RES«R0a.CF 

20  FORMAT  < * 0 S/A  <tE12»4»,0  R/A  'iE20.12,»  CORRELATION  PARAMETER 
1 ' • E2  0 . 1 2 I 

100  CONTINUE 
STOP 
END 
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APPENDIX  B 

METHODOLOGY  FOR  EXTRACTION  OF  EXPERIMENTAL  DAMPING  COEFFICIENTS 

In  the  case  of  a pendulum  such  as  the  one  shown  in  Fig.  B-l,  its  motion  is  the 
result  of  two  types  of  forces.  One  is  positional  in  nature  and  is  the  resultant  of  the 
gravitational  force  acting  on  the  pendulum  (at  its  center  of  gravity,  cg)  and  the  tension 
in  its  support  rod  (or  pendulum  arm).  As  shown  in  the  figures,  this  force  has  a magnitude 
of  Mg  sin  9 where  M is  the  mass  of  the  pendulum,  g is  the  acceleration  due  to  gravity, 
and  9 is  the  angular  position  of  the  pendulum  with  respect  to  the  Earth's  vertical.  In 
the  case  where  9 is  small,  sin  9 ~ 9,  and  the  restoring  force  is  seen  to  have  a magnitude 
of  Mgx/8  where  8 and  x are  apparent  in  Fig.  B-l. 

The  other  forces  acting  on  the  pendulum  are  frictional  in  nature  and  vary  with  the 
pendulum's  angular  velocity  (dfi/dt)  where  t is  time.  These  forces  arise  from  numerous 
sources.  The  dominant  ones  in  the  present  case  arise  from  air  resistance  and  magnetic 
forces  acting  at  the  center  of  magnetic  pressure.  Since  these  magnetic  forces  act  on  that 
portion  of  the  body  which  is  displaced  by  an  amount  (x’),  rather  than  x,  the  velocity 
(dx’/dt)  determines  the  magnetic  damping.  As  shown  in  Fig.  B-l,  dx'/dt  is  equal  to 
(Cm  /£)dx/dt.  where  8^  is  apparent  from  the  figure,  so  that  the  frictional  forces  must 
have  a magnitude  of  [(8m/8)Kmd  + CoKdx/dt)  where  Km  d is  the  magnetic  damping 
coefficient  and  Cp  is  predominantly  the  aerodynamic  drag  coefficient  but  includes  all 
other  effects  with  the  exception  of  magnetic  ones. 

The  equation  of  motion  of  the  pendulum  is  given  by  Newton's  second  law  so  that 

c“)-jr  (B-'> 

Since  the  pendulum  displays  damped  harmonic  motion,  the  solution  of  interest  is 
X = A|exp[-  (Kmd£m/8  + CD  )<t/2M)]j-  cos  (cut  - 5)  (B-2) 

where  A is  an  arbitrary  constant,  co  is  equal  to 

>Jg/8-  [(Kmd8m  /8  + CD  )/ 2 m]  2 
and  5 is  a possible  angular  displacement. 

It  is,  therefore,  possible  to  find  the  value  of  Km  d £\i  /£  + Cp  by  observing  the 
motion  of  the  pendulum.  The  damped  motion  is  typified  by  that  suggested  in  Fig.  B-2 


M 


dzx 

dtz 


Mgx  _ / K m D ! 
£ \ £ 
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so  that  it  is  possible  to  eliminate  concern  over  the  cosine  term  by  measuring  the 
displacement  (x)  at  times  when  cos  (a>t-5)has  the  value  of  unity;  i.e.,  at  the  extremes 
of  the  oscillatory  motion.  Then,  it  is  simply  a matter  of  measuring  the  magnitude  of 
x at  two  such  occasions  and  recording  the  time  interval  between  them  to  have  all 
information  needed  to  calculate  Kmd^m/2  + Cd- 

By  arbitrarily  assigning  the  zero  value  to  the  time  of  the  initial  measurement, 

X,  = A exp[-(KMD*WK +CD)t!/2M]|  ^ = A (B-3) 

the  weighing  factor  (A)  can  be  determined.  Then,, 

X2  = Xi  exp  [-  (Km  d /«  + CD ) At/2M]  (B-4) 


and  since  (Kmd&m/K  + Cd)  remains  the  only  unknown,  algebraic  manipulation  yields 


km  d 

e 


+ Cd  = ™ jir 

At  X! 


(B-5) 


The  coefficient  Cd  can  be  detennined  at  a time  when  Km  d is  zero  from  the  expression 


Cd  - 


2M 

At 


Xi 


(B-6) 
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NOMENCLATURE 

a Radius  of  conductive  shell  or  projectile  forebody,  m 

B Magnetic  induction,  Wb/m2 

b Distance  from  z-axis  to  guideway  line  currents,  m 

cm  Effective  center  of  magnetic  forces 

cg  Center  of  mass 

cp  Center  of  aerodynamic  pressure 

F Force  per  unit  length,  N/m 

fcoir  Correlation  parameter 

g Acceleration  of  gravity,  m/sec2 

h Length  of  conductive  shell,  m 

I Electric  current,  amp 

ICy  Magnetic  positional  force  constant,  N/m 

Km  d Magnetic  damping  force  constant,  N-sec/m 

C Moment  arm  for  aerodynamic  forces,  m 

Cm  Moment  arm  for  magnetic  forces,  m 

M Pendulum's  mass,  kg 

n nth  current  segment 

:R  Radial  position  of  projectile’s  axis  with  respect  to  guideway  axis,  m 

r Radial  distance  from  axis  of  projectile  or  other  reference  to  item  or  point  of 

interest,  m 
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S Shielding  ratio  (see  Eqs.  (1)  and  (2)) 

«.  Guideway  current  segment  length  (see  Fig.  1),  m 

t Time,  sec 

V Velocity  of  projectile's  cg,  m/sec 

W Energy  of  eddy  currents  induced  in  shell,  J 

x Cartesian  coordinate,  m 

y Cartesian  coordinate,  m 

Dispersion  of  projectile's  cg  from  z-axis,  m 

z Cartesian  and/or  polar  coordinate  (coincident  with  the  range  and  guideway  axis), 

m 

z Offset  z-axis,  m 

a Angle  of  attack,  rad 

P Angular  position  of  shell  referenced  to  x-axis,  rad 

7 Angle  between  V and  Z-axis,  rad 

S Phase  angle,  rad 

e Phase  lag  between  inducing  field  and  Field  of  eddy  currents,  rad 

M Permeability,  h/m 

p Volume  resistivity,  ohm-m 

C Area  resistivity,  ohm 

£2  Ohm 

co  Angular  frequency,  rad/sec 
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Subscripts 

D Damping  field,  currents,  forces,  etc. 

DRG  Drag 

g Gravity 

INF  Infinite 

M Magnetic  (positional) 

MD  Magnetic  damping 

M.H  Radial  M 

M/3  Angular  M 

MD  R Radial  MD 

M Dl?  Angular  MD 

n nth  segment 

o Refers  to  conditions  of  free  space 

P Positional  field,  currents,  forces,  etc. 

p Pressure 

Radial  component 

r Radial  direction  with  respect  to  currents 

s Source 

TRP  Transposed 

0 Angular  component 
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Superscripts 

Implies  an  offset  system  of  polar  coordinates 
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